Expectations for dihadron correlation measurements at the LHC 
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The suppression of high transverse momentum {Pt) inclusive hadron spectra in heavy-ion colli- 
sions as compared to the scaled expectation of high Pt hadron production in p-p collisions is the 
most direct manifestation of the interaction of hard partons with the soft bulk medium produced in 
heavy-ion collisions which is absent in p-p collisions. Yet the measured nuclear suppression factor 
Raa is a very averaged quantity and hence only a limited amount of information about the medium 
evolution and the nature of the interaction with the medium can be deduced from Raa. Measure- 
ments of hard back-to-back hadron correlations in 200 AGeV Au-Au collisions at RHIC have opened 
a new window to study the energy loss of partons in a medium in a more differential way and for a 
different distribution of in- medium pathlengths than in the case of Raa. In this work, we present 
an extrapolation of our results for back-to-back yields at RHIC energies to 5.5 TeV Pb-Pb collisions 
at the CERN LHC. We also discuss differences and similarities between the measurement at RHIC. 

PACS numbers: 25.75.-q,25.75.Gz 
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I. INTRODUCTION 

The suppression of single inclusive high Pt hadrons in 
heavy-ion colhsions as compared to the scaled expecta- 
tion from p-p collision has long been regarded as caused 
by energy loss of high pr partons into a dense partonic 
environment [ll, 0, S 0, a Q ■ ^^ RHIC, the nuclear sup- 
pression factor Raa for pions in central Au-Au collisions 
has been measured out to 20 GeV [2] and a factor ~ 5 
suppression observed. 

The hope of using hard probes such as the single hadron 
suppression in the context of heavy-ion collisions is to 
infer properties of the medium and its density evolu- 
tion from induced changes relative to the corresponding 
known hard process in p-p collisions. There is, however, 
growing evidence that single hadron suppression is not 
enough to unambiguously determine even the distribu- 
tion of energy/momentum shifts of partons traversing the 
medium |8|, |9|. However, if energy loss cannot be deter- 
mined reliably in momentum space, there is little hope 
to try to infer the QCD matter density distribution in 
position space. 

Back-to-back correlations of hard hadrons flO . Ill| are 
a more differential hard process. Due to the different 
geometrical averaging in the case of single hadron sup- 
pression vs. back-to-back suppression, one may hope to 
obtain information about the spatial distribution of dense 
medium from combining the two observables. While the- 
oretical studies for back-to-back correlations as measured 
at RHIC have been done [l2, llj, [iJ] , they seem to indi- 
cate that for RHIC kinematics the amount of additional 
information is not substantial, in essence various mod- 
els for the energy loss mechanism and for the density 
evolution which describe single hadron suppression also 
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perform well when compared with dihadron suppression. 
The notable exception is a class of energy loss models 
based on a purely linear dependence of the energy loss 
with pathlength — those seem to be strongly disfavoured 
by back-to-back correaltions [l3|- As suggested in \\^ . 
the reason why there is only little sensitivity to the QCD 
matter density distribution at RHIC kinematics may be 
that the lever-arm in momentum is not large enough to 
probe substantial shifts in parton momentum — for a 
steeply falling parton spectrum, even a moderate shift in 
parton momentum effectively resembles an absorption of 
partons, and this fact greatly reduces the sensitivity. At 
the LHC however where the partonic pT range is large 
this ceases to be a problem and consequently the suppres- 
sion of hard back-to-back correlated hadrons becomes a 
promising probe. 

In this paper, we aim to provide a baseline prediction for 
the per-trigger yield in hard back-to-back correlations. 
This complements a baseline prediction of the nuclear 
suppression factor Raa jl6| made within the same model 
framework. We compare with the calculation at RHIC 
kinematics and point out similarities and differences. 



II. THE FRAMEWORK 

As in [iJ, [13] we calculate the correlation strength of 
hadrons back to back with a hard trigger in a Monte- 
Carlo (MC) simulation. There are three important build- 
ing blocks to this computation: 1) the primary hard 
parton production, 2) the propagation of the partons 
through the medium and 3) the hadronization of the pri- 
mary partons. Only the step 2) probes medium proper- 
ties, and hence it is here that we must specify details for 
the evolution of the QCD medium and for the parton- 
medium interaction. Let us first discuss steps 1) and 3) 
which are common to the simulation in p-p and Pb-Pb 
collisions. 



A. Primary parton production 

In Ref. [13| it has been demonstrated that leading order 
(LO) perturbative Quantum Chromodynamics (pQCD) 
is rather successful in describing the PT-spectrum of in- 
clusive hadron production over a wide range in ^/s when 
supplemented with a -y/s-dependent iiT-factor to adjust 
the overall normalization. This factor parametrizes next- 
to-leading order effects. Since we are in the following 
only interested in ratios of PT-distributions, i.e. yields 
per triggered hadron, any factor independent of Pt drops 
out. Hence, in the following we use LO pQCD expressions 
without trying to adjust the absolute normalization. 
The production of two hard partons k, I with transverse 
momentum pT in LO pQCD is described by 



= 2_/'^i-^V^(^i''9^)^2/j7B(a;2,(9^)- 



da^i- 



i'j 



dt 



^ AB-^kl+X 

dp'^dyidy2 

"J 

(1) 

where A and B stand for the colliding objects (protons 
or nuclei) and yi(2) is the rapidity of parton k{l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction xi and a factorization scale Q ^ px is 
fi/Ai^ijQ^)- The distribution functions are different for 
the free protons [l^ ]l^ and nucleons in nuclei [20, |2l| . 
The fractional momenta of the colliding partons i, j are 



given by xi^: 



v^ 



(exp[±yi] -|-exp[±2/2])- 



Expressions for the pQCD subprocesses ^^—^ — (s, £, u) as 
a function of the parton Mandelstam variables s,t and 
ii can be found e.g. in [2^. By selecting pairs of k,l 
while summing over all allowed combinations of i, j, i.e. 
gg, gq, gq^qq,qq,qq where q stands for any of the quark 
flavours u, d, s we find the relative strength of different 
combinations of outgoing partons as a function oi pt- 
For the present investigation, we require yi = 2/2 = 0, i.e. 
we consider only back-to-back correlations detected at 
midrapidity. In the first step, we sample Eq. ([1]) summed 
over all fc, I to generate px for the event. In the second 
step we perform a MC sampling of the decomposition of 
Eq. ([1]) according to all possible combinations of outgoing 
partons fc, I at the px obtained in the first step. We thus 
end with a back-to-back parton pair with known parton 
types and flavours at transverse momentum px ■ 
To account for various effects, including higher order 
pQCD radiation, transverse motion of partons in the 
nucleon (nuclear) wave function and effectively also the 
fact that hadronization is not a collinear process, we fold 

,. ., . , AB-,kl + X 

mto the distribution J a ^ — :; — an intrinsic transverse 
momentum kx with a Gaussian distribution, thus cre- 
ating a momentum imbalance between the two partons 
as pxi + px2 = ^T- We use a value of kx = 2.7 GeV 
for RHIC kinematics [23[ and observe a sensitivity of less 
than 10% of the away side yield to changes of this number 
by ±50%. While kx is expected to rise at LHC condi- 
tions, in the momentum region we wish to explore, i.e. 



Px > 50 GeV there is no sensitivity of the calculation to 
the choice of intrinsic transverse momentum. 



B. Hadronization 

Before we can count hadrons above certain trigger or 
associate-momentum threshold, we have to convert the 
simulated partons into hadrons. More precisely, in order 
to determine whether there is a trigger hadron h above 
a given threshold, given a parton k with momentum px, 
we need to sample A'l~*'^{zi,px), i.e. the probability dis- 
tribution to find a hadron h from the parton k where h is 
the most energetic hadron of the shower and carries the 
momentum Px = Zi ■ px- 

In previous works |12l . [l3!| we have approximated this by 
the normalized fragmentation function Dk^h{z, n), sam- 
pled with a lower cutoff Zmin which is adjusted to the 
reference d-Au data. This procedure can be justified by 
noting that only one hadron with z > 0.5 can be pro- 
duced in a shower, thus above 2; = 0.5 the distributions 
Dk^h{z,Px) and Ai^'^{zi,px) are (up to the scale evo- 
lution) identical, and only in the region of low z where 
the fragmentation function describes the production of 
multiple hadrons do they differ significantly. However, 
for a fixed cut in the hadronic momentum Px, low z im- 
plies probing the high px part of a steeply falling parton 
spectrum, thus this region is suppressed. Empirically, 
sampling Dk—,h(z, Px) has provided a good description 
of the baseline [ij] ■ 

However, at the LHC energies, keeping the same lower 
cutoff Zmin is uot acceptable as it would imply that there 
is no perturbative production of hadrons in the range 
below, say, 8 GeV from a 100 GeV parton. Presumably, 
a viable solution would be to readjust the cutoff to a 
measured set of reference data from p-p collisions. In 
the absence of such a data set, we have chosen a different 
path by extracting A'l~^^ (zi , px) from shower simulations 
of quarks and gluons in PYTHIA [24] . 
This moreover has an additional advantage, since there 
is a subtle point with the scale dependence of the 
hadronization functions. While Dk^h{z,Px) is an em- 
pirical function where the scale is set by the produced 
hadron and the function is folded with the full parton 
spectrum (so that Px = zpx), A\~^^{zi,px) is an object 
where the scale is set by the given parton energy and the 
outcome of the folding is a distribution of hadron mo- 
menta. The latter is more appropriate for the situation 
realized in the MC simulation in which the partonic px 
is fixed before hadronization. 

In Fig. [T] we show the extracted A\~^^''{zi,px) compared 
with Dk^fjz, fj) from the KKP set of fragmentation 
functions [25[ in a probabilistic interpretation where we 
present the fragmentation function both at /i = Q.hpx 
and 11 — Px, allowing a direct comparison at zi = 0.5 
and zi = 1. 

In the region z < 0.5, the quantities Ai^'^(zi,px) and 
Dk^h{z, fj.) in Fig. [T] differ from each other as expected. 
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FIG. 1: (Color online) Comparison of the KKP [251 ] fragmen- 
tation function D{z,y?) for u-quarks and gluons into charged 
hadrons [K^ + h~) at two hadronic scales with the leading 
hadron momentum fraction probability density Ai{zi,pt) at 
partonic pr = 25 GeV as extracted from shower simulations 
in PYTHIA H. 



From the figure, it is also evident that in the region 
z > 0.5 where we expect agreement there remain differ- 
ences between the empirical fragmentation function and 
the shower development in PYTHIA which are beyond 
the uncertainties introduced by the QCD scale evolution. 
Especially in the case of the gluon fragmentation the dif- 
ferences are sizeable at higher z. We do not wish to enter 
a discussion of the physics origin of these differences at 
this point but rather accept Ai(zi,^) from PYTHIA as 
a trial ansatz which needs to be tested against data. 
Using Ai~^'^{zi,pt) on both outgoing partons from the 
primary pQCD vertex (see Eq. [T|) we find the two most 
energetic hadrons from both showers. Following the ex- 
perimental analysis procedure, out of the two we define 
the harder one to be the trigger hadron if it passes a 
set momentum cut, otherwise we ignore the event. The 
direction of the trigger hadron is taken to be the {—x) 
direction and the hemisphere corresponding to this direc- 
tion is referred to as 'near side', the opposite hemisphere 
is then defined to be the 'away side'. 
In order to compute the strength of hadrons correlated 
with a trigger on the away side, the leading contribution 
is found by considering the leading hadron produced from 
the away side parton for which Ai{zi,^) is sufficient. 
However, in order to determine the correlation strength 
on the near side, the production of (at least) one more 
energetic hadron in the shower needs to be considered 
as the leading contribution here corresponds to the trig- 
ger hadron itself. Thus, note that to the same level of 
approximation of the shower development the near side 
correlation is only found at next-to-leading (NL) frag- 
mentation whereas an away side correlation appears al- 
ready at leading fragmentation and receives a correction 
at next-to-leading fragmentation. 

In our framework, we treat this next-to-leading con- 
tribution through the conditional probability density 



A2{zl,Z2^^J) which for fixed zi corresponds to a prob- 
ability density, i.e. J dzidz2Ai{zi, fi)A2{zi,Z2, fJ.) — 1. 
This is conceptually similar to the dihadron fragmenta- 
tion function introduced in [26]. The whole hadron pro- 
duction in the shower arises in this language as a tower 
of conditional probability denities Aj^{zi, . . . ^Zn, ^) with 
the probability to produce n hadrons with momentum 
fractions zi , . . . z„ from a parton with energy /i being 

n?=iAj(zi,...Zi,^). 

In previous calculations [l^, llj] ^^ have modelled the 
next-to-leading conditional fragmentation probability us- 
ing the experimentally measured probability distribution 
Aj(zr) of associated hadron production in d-Au collisions 
[lOl [Il| as a function of Zt where Zt is the fraction of 
the trigger hadron momentum carried by the associated 
hadron. A factor 6{Ei — E'trigger — ^E — i^assoc) on the 
near side and 9{Ei — -Epunch — ^E — £'assoc) on the away 
side was included to enforce energy conservation. 
In this work, we use A2(zi,Z2,/z) as extracted from 
PYTHIA instead. By construction, Z2 < zi is satis- 
fied, and energy-momentum conservation (which the ex- 
tracted distribution fulfills automatically) implies Z2 < 
1 — zi, i.e. the region where A2(zi,Z2,/i) is nonzero in 
the (zi,Z2) plane is a triangular region from the edges 
zi = 0, zi = 1 where Z2 = to the center zi — 0.5 where 
z^"°^ = 0.5. 

The resulting distributions for u-quarks and gluons are 
shown in Fig. [51 



C. The medium model 

We describe the production and spacetime evolution 
of QCD matter in central Pb-|-Pb collisions within the 
framework of pQCD -I- saturation -I- hydrodynamics (the 
EKRT model) [23, 111, [2^. In this framework, primary 
parton production is computed by supplementing the 
coUinearly factorized pQCD parton production of Eq. (|T]) 
with the conjecture of saturation of produced gluons in- 
hibiting the growth of gluon number at transverse mo- 
menta below the saturation scale, pT < Psat- In cen- 
tral heaviest-ion collisions at RHIC and LHC, the com- 
putation remains in the perturbative regime as Psat ~ 
1.2 — 2.0 GeV ^ Aqcd- The average formation time of 
the QCD system, tq - l/^at = 0.17 - 0.2 fm/c at RHIC 
and LHC, and pQCD parton production in a central ra- 
pidity unit can then be computed, and thus the initial 
QGP densities estimated, see [2l, [2g| for details. 
With the assumptions of full thermalization and no ini- 
tial transverse flow at tq, we then apply relativistic, longi- 
tudinally boost-symmetric, azimuthally symmetric, lon- 
gitudinally and transversally expanding ideal hydrody- 
namics to describe the spacetime evolution of the formed 
QCD matter. We use the results obtained in [23], where 
the QGP is assumed to be an ideal gas of gluons and three 
flavours of massless quarks and antiquarks, while the 
hadron resonance gas consists of all hadronic states up 
to mass 2 GeV. A first-order phase transition takes place 




FIG. 2: (Color online) Conditional probability density A2{zi, Z2, f-i) for u-quarks (left) and gluons (right) for production of 
a next-to-leading hadron with momentum fraction Z2 given that a hadron with momentum fraction zi was produced in the 
shower at a parton energy scale /i = 25 GeV. 



at a critical temperature 165 MeV, and, with binary- 
collision initial energy-density profiles, the system de- 
couples at Tdoc = 150 MeV (see also the discussion of 
the decoupling dynamics and Tdcc in 31]). Illustrations 
of the phase boundaries, flow lines and transverse veloc- 
ity contours for central Au+Au coUisinons at RHIC and 
Pb+Pb at LHC, as well as extensive comparison with the 
hadron spectra measured at RHIC and further discussion 
can be found in |2g |. 

Within the framework we apply here, there are obvious 
theoretical uncertainties in the extrapolation to the LHC 
energies, which are related to the gluon shadowing (see 
p32 |) in the primary production, and to the QCD mat- 
ter equation of state and to the currently debated possi- 
ble dissipative effects [33, [SJ] in the matter evolution it- 
self. Accepting these uncertainties, however, the benefit 
in using the EKRT model is that being a closed frame- 
work it has predictive power: once Tdcc has been fixed at 
RHIC i2Q.], straightforward RHIC-tested predictions for 
the LHC hadron spectra can be made [2^, [SOl , and the 
evolution of QCD matter densities and flow velocities, 
the input for the present paper, obtained. 



D. Parton- medium interactions 

In the case of heavy- ion collisions, partons emerging from 
a hard vertex are surrounded by the soft bulk matter 
and may lose energy before hadronization. We compute 
the energy loss of a hard parton traversing a thcrmal- 
ized medium using the formalism by Baier, Dokshitzer, 
MuUer, Peigne and Schiff (BDMPS) which assumes 
that energy is predominantly lost due to medium-induced 
radiation. This formalism is cast into the form of en- 
ergy loss probability distributions, so-called quenching 
weights in 36\. This is a convenient formulation for the 



purpose of a MC simulation. 

The key quantity characterizing the energy loss induced 
by a medium with energy density e in the BDMPS formal- 
ism [2] is the local transport coefficient q{T,r]s,r) which 
characterizes the squared average momentum transfer 
from the medium to the hard parton per unit path- 
length. Since we consider a time-dependent inhomoge- 
neous medium, this quantity depends on proper time 
T = y/t^ — z^, spacetime rapidity r/s — ^\n |^ , cylin- 
drical radius r and in principle also on azimuthal angle 
(f), but for the time being we focus on central collisions 
only. 

The transport coefficient is related to the energy density 
of the medium e as 



q{T,r]s,r) = i4r-2- [e(T, 77^, r)]'^/^(coshp— sinhpcosa) (2) 

with K = 1 for an ideal quark-gluon plasma (QGP) [33| 
Here, p is the local flow rapidity with angle a between 
flow and parton trajectory [38, 39]. The parameters e 
and p characterize the medium and must be supplied by 
the underlying QCD medium evolution model. 
In the following, motivated by Ref. [33| we assume the 
proportionality constant K to remain unaltered in differ- 
ent phases of the medium. We treat K as an adjustable 
parameter (for a detailed discussion see [ij]) which how- 
ever, once fixed, is independent on y/s of the heavy-ion 
collision and assumes the same value at RHIC and LHC 
energies while the medium density given by e(r, 773,r) 
changes substantially. 

Given the local transport coefficient at each spacetime 
point, a parton's energy loss depends on the position of 
the hard vertex at ro = {xo, yo) in the transverse plane at 
T = and the angular orientation of its trajectory (jj (i.e. 
its path through the medium). In order to determine the 
probability P(A_E, E)path for a hard parton with energy 



E to lose the energy AE while traversing the medium 
on its trajectory, we make use of a scaling law [4Qj] which 
allows to relate the dynamical scenario to a static equiva- 
lent one by calculating the following quantities averaged 
over the jet trajectory ^(r) : 



as a function of the jet production vertex Tq and its angu- 
lar orientation (/). We set q = whenever the decoupling 
temperature of the medium T = Tp is reached. 



and 



Jo 



{qL){vo,<p)^ d^m 

Jo 



(3) 



(4) 



Using the numerical results of [3g|, we obtain P{AE)path 
for ujc and R = 2uj'^/{qL) as a function of jet production 
vertex and the angle corresponding to 



P{AE) 



path 



E: 

n=0 



n „ 



dI{oJi) 
dio 



s(ae-J2^ 



exp 



duj 



dl_ 

dio 



(5) 



which makes use of the distribution w^ of gluons emit- 
ted into the jet cone. The explicit expression of this 
quantity for the case of multiple soft scattering can be 
found in [3a|. Note that the formalism of 36] is defined 
for the limit of asymptotic parton energy, hence the prob- 
ability distribution obtained P{AE)path is independent 
oi E. 

The initially produced hard parton spectrum and, con- 
sequently, the number of hard vertices in the (x, y) plane 
(where the 2:-axis is given by the beam direction) in an 
A — A collision at fixed impact parameter b, are propor- 
tional to the nuclear overlap. 






->/+x 



dp^dyf 



dp^dyj 



(6) 



where TaaC^) is the standard nuclear overlap func- 
tion. We define a normalized geometrical distribution 
P{xo,yo) for central collisions as 



PixoiVo) 



Taa{0) 



(7) 



where the thickness function is given in terms of the nu- 
clear density yOA(r, z) as TA{r) = / dzpA{Y, z). 
For a given event in the MC simulation, we first sample 
Eq. ([7]) to determine the vertex position in the transverse 
plane. From this point, we propagate the parton through 
the medium and we evaluate Eqs. pi4p along the path, 
on which the medium properties enter via Eq. ^. Com- 
puting the corresponding P{AE)path for each parton, we 
sample this distribution to determine the actual energy 
loss for near and away side parton in the event. In the 
case AE > Eparton we consider the parton to be ab- 
sorbed by the medium and not contributing to high px 
yield. 



III. RESULTS 

A. Comparison with 200 AGeV d-Au data 

The first step is to test the MC simulation without a 
medium to see how well LO pQCD combined with the 
first two terms in the hadronization probability distribu- 
tion expansion, i.e. Ai(zi,^) and ^2(21,22,/^) describes 
the data. For this purpose, we compare in Fig. [3] with 
the correlation strength measured in d-Au collision by 
the STAR collaboration [10, llll. The data have been 
taken for a trigger hadron between 8 and 15 GeV and in 
two associate yield bins (on both near and away side) of 
4-6 and 6-1- GeV. In order to simulate d-Au collisions, we 
insert a nuclear parton distribution function into Eq. ([T]) 
and assume that no medium is formed and that nuclear 
matter is not dense enough to induce significant energy 
loss. 

As apparent from the figure, the calculation is able to 
describe the data generally rather well, however in the 
4-6 GeV momentum bin on the away side it overshoots 
the measured data. We have re-run the simulation tak- 
ing into account only the leading fragmentation term 
Ai (zi , fi) and find that this contribution alone seems to 
do well with the data. 

Since the trigger is always chosen to be the most en- 
ergetic hadron in the simulation, there is considerable 
bias towards the trigger parton being a quark, as the 
quark fragmentation is harder (cf. Fig. [1]). In the kine- 
matical range probed, the q{q)g -^ q{q)9 subprocess is 
rather important, thus there is a bias towards a gluon 
being the away side parton. We note from Fig. [1] that 
the discrepancy between the gluon fragmentation func- 
tion and Af~^ (^IjM) is larger than for quarks. Quite 
likely, this discrepancy also propagates into the extrac- 
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FIG. 3: (Color online) Yield per trigger-hadron between 8 and 15 GeV in d-Au collisions at 200 AGeV on near side (left panel) 
and away side (right panel) as obtained by the STAR collaboration [lOl . [ll| and compared with the MC simulation described 
in this work for two associate momentum bins of 4-6 GeV and 6-1- GeV. The different data points have been shifted artificially 
along the i-axis for clarity of presentation. The right panel also includes the result neglecting the NL fragmentation of the 
away-side parton. 



tion of v4f ^ (zi , Z2 7 m) from PYTHIA and is the reason 
that the present calculation shows differences here rela- 
tive to the data (and to previous results [12, [l3| ) • 
The eventual solution to this discrepancy would be to 
find a different way to extract Ai{zi, ^) and A2{zi, Z2, /i) 
than PYTHIA shower simulations. However, since the 
overall description of the data seems to be reasonable, we 
proceed with the hadronization procedure as introduced 
above, but assign an uncertainty to the simulations of 
the heavy-ion case which we estimate by calculating the 
results with and without ^2(2:1, 22, m) on the away side. 



produces the right order of magnitude to agree with the 
data. A different possibility is a sizeable contribution of 
coalescence processes [4l|, |42] in the Au-Au case to the 
near side hadron yield which is not included in the cal- 
culation. If coalescence processes are relevant, one would 
expect to systematically underestimate the data. Given 
the large uncertainty of the baseline d-Au measurement 
and the theoretical uncertainty in the shower modelling 
using PYTHIA, it is unfortunately impossible to make 
a definite statement what the relevant physics process 
underlying the discrepancy is. 



B. Comparison with 200 AGeV central Au-Au data 

In the next step, we test the MC simulation including 
in-medium energy loss with a soft medium background 
evolution as provided by a hydrodynamical calculation 
for 200 AGeV central Au-Au collisions. The parameter 
K in Eq. ^ is fixed by the single hadron suppression for 
this system [13|, thus the computation proceeds without 
any further free parameters. The resulting yields per 
trigger as a function of associate hadron momentum for 
near and away side are shown in Fig. [H 
We observe that the yield on the near side is almost iden- 
tical to the calculated yield in d-Au collisions and system- 
atically below the data. One way of resolving the discrep- 
ancy while retaining agreement with the d-Au data could 
be the inclusion of NNL fragmentation processes. Due to 
the large errors in the d-Au measurement, this would 
not destroy the agreement with the baseline. While a 
full computation of the NNL conditional fragmentation 
probability distribution ^3(2:1, 2:2, Z3, //) is very involved, 
we have checked that the approximation 



As,(zi,Z2,Z^,y) « ^2(2:1 -1-22,23,^)^(^2 - 23) (8) 



On the away side, the calculation shows good agreement 
in both the 4-6 GeV and the 6-|- GeV momentum bin. 
However, given that the calculation overshoots the d-Au 
baseline in the lower momentum bin, one has to conclude 
that the agreement in this bin is accidental. 



This is clearly confirmed by omitting NL fragmentation 
on the away side (which restores agreement with the base- 
line): doing so results in an underprediction of this bin 
just as observed in our previous calculation \'S\ which 
is consistent with the observation that coalescence pro- 
cesses [4l|, I43, |4j1 play a substantial role for hadronization 
in this momentum region. 



In Fig. [3 we show the distribution of zt = 
passoc ^ p^igger ^^ ^j^^ away side in d-Au and central 

Au-Au collisions and compare with data. We observe 
a general tendency that the calculation slightly underes- 
timates the data in the high zt region. This is consistent 
with the observation that there is a trend that the cal- 
culation may underestimate the higher associate momen- 
tum window slightly. However, within errors, the results 
are compatible. 
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FIG. 4: (Color online) Yield per trigger hadron between 8 and 15 GeV in central Au-Au collisions at 200 AGeV on near side 
(left panel) and away side (right panel) as obtained by the STAR collaboration [lOl . [lll | and compared with the MC simulation 
with a hydrodynamical evolution for the soft medium background described in this work for two associate momentum bins of 
4-6 GeV and 6+ GeV. The different data points have again been shifted artificially along the a;-axis for clarity. The right panel 
also includes the result neglecting the NL fragmentation of the away side parton. 
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FIG. 5: (Color online) zt distribution of away side associated 
hadrons measured in 200 AGeV d-Au and central Au-Au col- 
lisions in comparison with the MC simulation results. Here, 
'baseline calculation' (see text) denotes the simulation with- 
out a medium and 'hydrodynamics' the simulation carried out 
using a hydrodynamical evolution model for the QCD matter 
produced. 



C. Prediction for 5.5 ATeV central Pb-Pb collisions 

We compute the expected correlation strength pattern 
for an energy of 5.5 ATeV for p-p coUisions and Pb-Pb 
colhsions as expected to take place at the CERN LHC. 
For the trigger range, we assume 50-70 GeV to probe a 
regime which can be accessed at LHC but is beyond the 
reach of the RHIC experiment. 

Since our formalism includes shower development only 
to the second term A2 (zi , 2:2 , /i) , the associate hadron 
momentum cuts cannot be too low, otherwise the main 
contribution to hadron production cannot be assumed to 
be dominated by the leading and next to leading hadron. 
Thus, we start considering associate hadron production 
for 20 GeV and above up to the trigger energy in 5 GeV 



momentum windows. In order to account for the uncer- 
tainty with regard to the NL fragmentation on the away 
side as discussed earlier, we carry out the simulation both 
including and excluding this contribution to provide an 
error estimate. 

In Fig. [S] we show the resulting yield per trigger on the 
near and away side. Let us first focus on the away side 
(right panel). Regardless if one uses NL fragmentation 
contributions or not, there is a suppression of the away- 
side yield when going from p-p to Pb-Pb collisions of 
about a factor 3. This is comparable with the factor 4 
suppression observed at RHIC. The similarity is not quite 
unexpected, since in the range of the trigger momentum 
the same model predicts a single hadron suppression fac- 
tor of about the magnitude seen at RHIC ,16|. The un- 
derlying reason is that the ratio between the typical scale 
of energy loss vs. the momentum scale probed is similar 
and the change in characteristic length scales between 
central Au-Au and central Pb-Pb is negligible. Up to a 
momentum scale of about 35 GeV, there is an uncertainty 
associated with NL fragmentation, for higher momenta 
the uncertainty goes away . 

Let us now focus on the near side (left panel in Fig. [H]): 
Here we observe that the yield per trigger is enhanced 
by 10-20% in Pb-Pb relative to to p-p collisions. At 
first glance, this seems to be surprising. However, one 
should keep in mind that the model still includes a strong 
suppression of trigger hadrons, only the per-trigger yield 
is enhanced. These findings can be reconciled by the 
observation that the dense medium at LHC acts like a 
gluon filter and in essence enhances the relative abun- 
dance of quarks escaping the medium to gluons escaping 
the medium. The underlying reason is that gluons, due to 
their different color charge, couple more strongly to the 
medium than quarks and consequently lose more energy. 
However, due to the harder leading and NL fragmen- 
tation, quark jets on average produce harder associated 
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FIG. 6: Yield per trigger hadron between 50 and 70 GeV in p-p and central Pb-Pb collisions at 5.5 ATeV on near side (left 
panel) and away side (right panel) as obtained in our MC simulation with a hydrodynamical QCD medium evolution as a 
function of associate momentum bin. The right panel also includes the result neglecting the NL fragmentation of the away-side 
parton. 
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FIG. 7: (Color online) Yield per trigger hadron between 50 
and 70 GeV in p-p and central Pb-Pb collisions at 5.5 ATeV 
on the away side as obtained in our MC simulation. The yield 
in both Pb-Pb cases is scaled with a constant to agree with 
the first point of the p-p computation. The distortion of the 
momentum spectrum of recoil hadrons due to energy loss is 
clearly visible. 



yield than gluon jets, which is precisely what is seen in 
the simulation. 

Finally, in Fig. [7] we take a closer look at the niomentum 
spectrum of hadrons recoiling from the trigger. To first 
order, one finds the reduction by a factor three mentioned 
above when going from p-p to Pb-Pb collisions. This 
would be consistent with complete absorption of 2/3 of 
all partons whereas 1/3 passes the medium without any 
energy loss. However, as a closer comparison reveals, 
there is a momentum dependent distortion of the spectra 
in heavy-ion collisions which is associated with events in 
which an away-side parton emerges from the medium af- 
ter a finite energy loss AE. From this pr dependent shift 
of the spectrum, one can in principle derive constraints 
about the averaged energy loss probability P{AE, E) and 
hence about the medium density distribution. 



IV. DISCUSSION 

We have made a baseline prediction for the expected per- 
trigger yield of associate hadrons as a function of trans- 
verse momentum in central Pb-Pb collisions at 5.5 ATeV 
for a trigger energy of 50 - 70 GeV. We have demon- 
strated that the effect of the medium is expected to leave 
a characteristic imprint in the observables. 
On the near side, we expect enhancement of the associate 
hadron yield as compared to the p-p case due to the 
fact that energy loss of gluons is stronger, hence there is 
some bias towards quarks leading to a trigger hadron, and 
the fact that the NL fragmentation for quarks is harder 
leads to a slight enhancement in the relevant momentum 
region. 

On the away side, we expect a substantial suppression of 
the yield as compared to the p-p case, although the sup- 
pression we find is in relative terms smaller than what 
is currently observed at RHIC. We have further shown 
that the shape of the momentum spectrum of away-side 
hadrons is not equal to the momentum spectrum ex- 
pected for p-p collisions. This is because the medium 
does not only act by absorbing partons, but one gets a 
chance to observe the characteristic shift in parton energy 
induced by the medium. This allows a more detailed view 
on the geometry-averaged energy loss probability distri- 
bution. In this respect, the observable shows a potential 
which may not be accessible at RHIC kinematics. 

There are multiple uncertainties associated with this pre- 
diction at each step of the MC simulation. In the LO 
pQCD calculation, both NLO corrections and the choice 
of the nuclear parton distribution functions result in 
an uncertainty, especially as the role of gluon shadow- 
ing cannot be reliably extrapolated to the LHC regime. 
However^ we have tested that the difference between the 
NPDF ,203 and the EKS98 [2l| nuclear parton distribu- 
tion sets is smaller than the statistical errors of the MC 
simulation. While the yield of partons at given pT is sen- 



sitive to the nuclear parton distribution, the associated 
per-trigger yield is not since many uncertainties cancel 
in the ratio. 

There is likewise an uncertainty associated with the 
choice of the hadronization scheme. We have tried to 
show the magnitude of this uncertainty by comparing 
with available RHIC data and by indicating the relative 
magnitude of the leading and next-to-leading fragmenta- 
tion term in the prediction. 

A large uncertainty must also be assigned to the hydrody- 
namical prediction of the underlying medium evolution. 
At present, neither entropy production and the distribu- 
tion in the initial state nor the equilibration process or 
the equation of state governing a hydrodynamical expan- 
sion are precisely known. Yet all these quantities would 



leave a characteristic signature on the geometrical bias 
induced by the requirement to obtain a trigger. 
Keeping in mind these uncertainties, the prediction pre- 
sented here should in any case show the right order of 
magnitude of the correlation strength. Furthermore, we 
believe the essential effects are described qualitatively 
correct, given that the radiative energy loss picture is 
valid. 
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